﻿ Lucrarea nr 5 Modelarea matematică a sistemelor de transmisiuni în banda de bază fără interferenţă intersimbol Scopul lucrării: se urmăreşte familiarizarea studenţilor cu problematica modelării matematice a unor fenomene şi sisteme fizice Se are în vedere o modelare pe calculator a unui sistem de transmisiuni de date în bandă de bază, modelare ce presupune o prelucrare numerică, (în timp discret), în timp ce “modelul” matematic stabilit este o funcţie de variabile continue 1 Modelarea sistemelor de transmisiuni în banda de bază Teorema lui Nyquist Prin transmisiune în banda de bază înţelegem acea transmisiune care utilizează porţiunea din apropierea frecvenţei zero din banda de frecvenţă a canalului Semnalele în banda de bază, impulsuri dreptunghiulare, au un spectru de frecvenţe care ocupă, teoretic, o bandă infinită Însă banda de frecvenţe utilizabilă la o transmisie de date este limitată atât datorită canalului cât şi din considerente econommice În figura 1 se prezintă schema bloc a unui sistem de transmisiuni în banda de bază: z(t) a(t) s(t) s’(t) r(t) {âk} {ak} antionare/ Filtru de Filtru de Eş MIA Canal de com ie Decizie Emisie Recepţ ) G(ω C(ω) E) GR(ω tk )=GE(ω) C(ω) GR(ω) G(ω Figura 1 Schema bloc a unui sistem de transmisiuni în banda de bază Informaţia ce trebuie transmisă este conţinută în succesiunea de impulsuri a modulate în amplitudine, a(t), de secvenţa de date de la intrarea sistemului { } k ) dL, ,d,d13−±±±, L fiind 2 pentru semnale Amplitudinile posibile sunt ( binare La emisie şi recepţie se utilizează filtrele cu funcţiile de transfer ( Gω ) E şi respectiv () Gω Semnalul recepţionat r(t) este prelucrat de către receptor R pentru a se stabili amplitudinea sa şi, în consecinţă, simbolul ce a fost emis 1 C ω este funcţia de transfer a canalului ideal Răspunsul în frecvenţă al () întregului lanţ de transmisiuni (emiţător-canal-receptor), ( Gω, este: ) GG CGωωωω=⋅⋅ (1 1) () () ()() ER numit şi răspunsul în frecvenţă al canalului echivalent Deorece banda de frecvenţe utilizabilă la o transmisie de date este limitată, rezultă că acest răspuns în frecvenţă al canalului echivalent trebuie să fie o funcţie de bandă limitată, ce impune însă şi o limită maximă a vitezei de transfer a datelor Teorema lui Nyquist afirmă că într-un canal echivalent cu un filtru trece-jos a de tăiere F, este posibilă transmiterea semnalelor binare ideal, cu frecvenţ independente, cu o viteză de semnalizare mai mică decât 2F simboluri/secundă , fără interferenţă între simboluri Figura 2 a prezintă funcţia de transfer de tip “filtru trece-jos ideal” (cu fază nulă), () G ω, în timp ce figura 2 b reprezintă răspunsul la impuls, (gt, ) corespunzător a) b) Figura 2 a) Răspunsul în frecvenţă al filtrului trece-jos ideal; b) Răspunsul la impuls corespunzător Expresiile analitice ale celor două funcţii sunt: 1 ⎧ Ftsin2π Fωω π≤= 2( ⎪ 0) Gω=, ()gtct==sinω (1 2) () F⎨ 2() 0 Ft2π ωω> 0 ⎪ 0 ⎩ Considerând o succesiune de impulsuri pe linie, decalate în timp cu durata T, condiţia de interferenţă nulă între simboluri se poate scrie sub forma: *1 ,ℵ ∈∀⋅=kkT (1 3) 2F 2 1 din care rezultă viteza de transmitere maximă : vF ==2 max T min Considerând un mesaj a(t), ce constă din N simboluri consecutive, ce apar la intervale de timp T, acesta are expresia: − 1N () −δ=nTta)t(a n (1 4)∑ =0n Semnalul recepţionat, r(t), în cazul unui canal ce nu introduce zgomot peste semnalul util (z(t)=0), are expresia: 1 N− rt a gt nT()=⋅− (1 5) () n∑ 0 n= Condiţia de interferenţă nulă intersimbol presupune ca eşantionul rrkT= k(), să fie dependent doar de simbolul a ş i de simbolurile ank, ≠ Notând ki nu şn ggkT=, eşantionele funcţiei răspuns la impuls g(t), condiţia de mai sus se () k verifică pentru: ≠0, (1 6) gnk=∀ nk− Se observă, analizând relaţia (1 2), că în condiţiile din relaţia (1 3), se verifică întotdeauna condiţia de interferenţă nulă intersimbol din ec (1 6) 2 Simularea pe calculator a trecerii unei secvenţe de date printr-un canal echivalent de tip “trece-jos ideal” Pentru a putea realiza o simulare a transmisiei cu ajutorul calculatorului, plecând de la filtrul global de tip “trece-jos ideal”, trebuie să găsim un model matematic discretizat (în timp discret), echivalent cu modelul în timp continuu stabilit În acest scop vom utiliza metoda invarianţei răspunsului la impuls 2 1 Metoda de echivalare a sistemelor în timp continuu cu sisteme în timp discret prin invarianţa răspunsului la impuls Metoda presupune considerarea ca şi răspuns la impuls al sistemului gn gnT=, obţinută prin eşantionarea discret echivalent secvenţa [ () ] de răspunsului la impuls g(t) al sistemului în timp continuu, cu un pas de teorema eşantionării În aceste condiţii Transformata eşantionare Te, ce verifică 3 Fourier în Timp Discret, (TFTD), a secvenţei [gn, (GΩ, se exprimă în ] d)d Gω conform: funcţie de răspunsul în frecvenţă ( ) ⎞⎛ 21 π ⎟⎜ () −ω=Ω ∑ dkGG (2 1) ⎟⎜ TTT ⋅ω=Ω e ⎠⎝ ∈Zkee , cu perioada π2 Ea este o funcţie periodică în Ω spunsul unui sistem în timp discret la o secvenţă de intrare []nx 2 2 Ră Notând cu [] ny, secvenţa de la ieşirea unui sistem în timp discret, cu δn , dat de secvenţa [], la o secvenţă de răspunsul la impulsul unitate, [ ng ] d intrare [] nx, dacă sistemul este liniar şi invariant în timp, avem: ∞ [] [ ] −⋅=kngkxny (2 2) [] ∑ −∞= k ceea ce reprezintă o sumă infinită de termeni Însă, în general, avem de-a face cu semnale de durată finită, caz în care indicele k din ecuaţia (2 2) ia valori dintr-un interval finit ] ny este oferită de O cale mai rapidă de calcul a valorilor secvenţei [ utilizarea algoritmului Transformatei Fourier Rapide (notat TFR sau FFT-Fast Fourier Transform), algoritm ce presupune un număr minim de operaţii necesare De fapt, aplicâd TFTD, ecuaţiei (2 2) rezultă: ()() (2 3) () Ω⋅Ω=ΩXGY d ][] ΩY şi()ΩX reprezintă TFTD a secvenţelor [ny şi respectiv nx unde () 2 ) cu kπ⋅=Ω Algoritmul TFR permite calculul rapid a valorilor ( kYΩk, M unde M este de obicei o putere a lui 2 Calculul rapid a acestor valori se face conform: () ()() ⋅Ω=Ω, (2 4) kkdkXGYΩ ] ng şi [nx deci utilzând eşantioane ale spectrelor secvenţelor [ d] Aplicând Transformarea Fourier Rapidă Inversă, (TFRI sau IFFT), secvenţei de ] se obţine secvenţa răspuns, [ny valori spectrale () kYΩ 4 2 3 Generarea cu calculatorul a unei secvenţe de date aleatoare Avem în vedere simularea comportării unui sistem în timp continuu, caracterizat prin funcţia răspuns în frecvenţă şi răspunsul la impuls date de ec (2) Această analiză presupune o secvenţă de date de intrare de forma dată de ec (4) Ne vom rezuma doar la o transmisie binară (deci cu doar două nivele posibile ale semnalului de date), unipolară, (nivele normalizate 0 şi 1 pentru 0 şi respectiv 1 logic), şi bipolară, (nivele normalizate -1 şi 1 pentru 0 şi respectiv 1 logic) În acest scop vom găsi sistemul discret echivalent al sistemului considerat, generând şi o versiune discretă a semnalului de intrare a(t) Am văzut, §2 1, că echivalarea sistemelor are la bază o eşantionare cu o se aplică şi perioadă T a unui semnal în timp continuu O procedură e similară semnalelor de intrare/ieşire Astfel, căutăm semnalul în timp discret []() , semnal ce se aduce la intrare sistemului discret echivalent ednTana= Ţinând cont de structura particulară a semnalului de intrare a(t), şi anume N simboluri binare de durată T, alegem perioada T astfel încât să e nu rezulte “pierdere sau denaturare de informaţie” (scăpare de bit sau sondarea diferiţilor biţi în număr diferit de puncte) Acest lucru se asigură alegând raportul TT L= eb un număr întreg, de obicei o putere întreagă a lui 2 Considerând un număr N de biţi, tot o putere întreagă a lui 2 rezultă o lungime totală a secvenţei numerice, MNL=⋅, o putere întreagă b a lui 2 Figura 3 O secvenţă de impulsuri unipolare ] na în care apare o valoare de bit la Se obţine astfel o secvenţă numerică [ d intervale de timp egale cu L b Dacă această valoare (0 sau 1 in cazul unipolar şi -1 sau 1 în cazul bipolar) se generează în mod aleator, (prin pe utilizarea unori funcţii ce generează valori aleatoare), obţinem o secvenţă de date de N simboluri succesive aleatoare Un exemplu de o astfel de secvenţă se prezintă în figura 3, unde s-a considerat un 5 4 număr de N=8 simboluri binare unipolare, cu 162=, rezultând M=128, =L b ] na) (lungimea totală a secvenţei numerice [ d 2 4 Răspunsul în frecvenţă al sistemului discret echivalent tip “filtru trece-jos ideal” Acesta are la bază ecuaţia (1 7), rezultând o funcţie periodică, cu perioada 2⋅π în variabila Ω Pulsaţia de tăiere a filtrului în timp continuu fiind antionare impune verificarea condiţiei ππTF=⋅2(vezi fig 2 a), teorema de eş l TT1≥ Astfel, se justifică o alegere de forma TT L l== ≥20, Expresia eeb răspunsului în frecvenţă al sistemului discret echivalent este: π ⎧ LΩ≤ b L⎪ b⎪ π2, (2 5) GGGΩΩΩ==+ ()() ( ) ddd⎨ π π0 Ω ⋅+ ωω α () 001 ⎪ ⎩ unde α este un parametru variabil între 0 şi 1, denumit şi exces de bandă Funcţia răspuns la impuls corespunzătoare, g(t), este: ttπ απ sin cos ∞ TTjt1ω (3 2) gt G e d==⋅ ωω ()() ∫ π t22π t⎛ 2α ⎞ −∞ 1 − ⎜⎟ T ⎝⎠ T ce verifică condiţia de interferenţă nulă intersimbol: gn kT n k−=∀≠0, (3 3) () () Gω pentru α=05 În figura 1 se prezintă funcţia răspuns în frecvenţă ( ) Figura 1 Răspunsul în frecvenţă al unui filtru de tip “cosinus-ridicat” Alegând α =1 se obţine funcţia de transfer global ce poartă numele de “filtru de tip cosinus pătrat” Relaţia (3 2) devine: ππω ⎧ 2 ωω≤ ⎪ 02cos (3 5) G ω = ωω () ⎨ 004 ⎪ ωω> 002 ⎩ Funcţia răspuns la impuls are expresia: π ttπ sin cos ∞ TTjt1ω (3 6) gt G e d==⋅ ωω ()() ∫ t22π 2π t⎛ ⎞ −∞ 1 − ⎟ ⎜ T ⎠ ⎝ T 8 ce verifică deasemenea condiţia de interferenţă nulă intersimbol, ec (3 4) Figura 2 a prezintă funcţia răspuns în frecvenţă, ( Gω (ec (3 5)), în timp ce ) gt (ec (3 6)) figura 2 b prezintă funţia răspuns la impuls corespunzător, ( ) a) b) Figura 2 a) Răspunsul în frecvenţă al filtrului “cosinus-pătrat”; b) Răspunsul la impuls corespunzător 3 1 Procedura de simulare a lanţului de transmisiuni în banda de bază Procedura de simulare implică: 3 1 1 Discretizarea semnalului de intrare (at, prin eşantionare cu pasul ) T Se obţine astfel o secvenţă numerică [an anT=, unde se alege () e]de l TT L l== >20,,, rezultând o lungime totală MNL=⋅ a secvenţei ebb Alegând numărul de simboluri de transmis N ca fiind o putere întreagă a lui 2, va rezulta o lungime totală a secvenţei, M, tot o putere întreagă a lui 2 3 1 2 Echivalarea sistemului în timp continuu cu un sistem în timp discret pe baza metodei de echivalare prin invarianţa răspunsului la impuls Se găseşte răspunsul în frecvenţă al sistemului numeric echivalent, ( GΩ, funcţie ce, în ) d l condiţiile TT L l== 20, (ce respectă teorema eşantionării), va avea eb > expresia: LΩ 22π b⎧ LΩΩ cos = ⎪ 0 b⋅≤ GLGGΩ =+,42π = (3 7) () () ( ) ⎨ dbddΩΩ ⎪ ΩΩ ω ⎩ T Funcţia răspuns la impuls, g(t), corespunzătoare, verifică condiţia de )() kn,Tkng≠∀=−0 interferenţă nulă intersimbol ( Considerând un mesaj a(t), ce constă din N simboluri consecutive, ce apar la intervale de timp T, acesta are expresia: − 1N () −δ=nTta)t(a n (4 4)∑ 0n = tg ră irea acestuia se Notând cu () Espunsul la impuls al filtrului de emisie, la ieş obţine semnalul: −1N −⋅=nTtga)t(s (4 5) () ∑ En 0n = semnal ce urmează a fi transmis prin canal Canalul fiind idealizat, deci nu atenuează semnalul, nu-l întărzie sau dispersează, şi ţine cont doar de zgomotele ce le introduce, la intrarea filtrului de recepţie vom avea semnalul: () () tztst's+= (4 6) () relaţie ce pune în evidenţă cele două componente ale unui semnal ce trece printr- un canal zgomotos: componenta de semnal util, s(t), şi componenta de zgomot, z(t) Dacă filtrul de recepţie este un sistem liniar atunci şi semnalul de la ieşirea acestuia va fi compus din două componente, una de semnal util şi cea de zgomot Astfel, semnalul recepţionat, r(t), se poate scrie: () () tztrtr+= (4 7) () u0 tr reprezintă spunsul filtrului de unde () u componenta de semnal util, deci ră ) tz este componenta de zgomot, deci răspunsul recepţie la semnalul s(t), iar ( 0 12 filtrului de recepţie la zgomotul z(t) din canal Componenta de semnal util are −1N expresia: () −⋅=nTtga)t(r (4 8) ∑ nu 0n = Semnalul recepţionat r(t) este eşantionat, iar pe baza eşantionului rrkT= k() se ia decizia asupra simbolului binar a transmis Semnalul recepţ kionat, r(t), este eşantionat la momentele de timp tkT= antioanelor obţinute k, iar pe baza eş = rrt=, se decide asupra simbolului transmis a Cum (), () kkkkuukakTrr= decizia corectă la recepţie presupune ca eşantionul de zgomot () kTzz=, în k00 modul, să nu depăşească un anumit prag (1/2 în cazul unei transmisii unipolare şi 0 în cazul unei transmisii bipolare) Se defineşte raportul semnal/zgomot, RSZ, la ieşirea filtrului de recepţie la momentul t ca fiind: 0 2 () 0utr RSZ = 0 (4 9) P 0z unde P este puterea zgomotului de la ieş irea filtrului de recepţie 0z 2 În cazul unui zgomot alb gaussian din canal, de putere σ i, puterea zgomotului la ieşirea filtrului de recepţie este: ∞ 2221 ωω⋅σ=σ=dGP (4 10) () Riz0 ∫ 0 2 π ∞− Semnalul util recepţionat, la momentul 0=t, se poate scrie: 0 ∞ 1 () () ωω⋅ω=dGSr0 (4 11) () Ru ∫ π 2 ∞− ) ωG, care să maximizeze raportul Se caută filtrul de recepţie optim, ( R semnal/zgomot la ieşirea filtrului de recepţie la momentul 0= 0t Acest filtru se găseşte ca fiind unul cu expresia răspunsului în frecvenţă: * () () ω⋅=ωSkG (4 12) R deci, filtrul de recepţie este adaptat la forma semnalului din canal Necunoscând însă, apriori, informaţia ce se va transmite prin canal, cerinţa (4 12) este dificil de realizat O soluţie mai puţin pretenţiosă, insă doar suboptimală, acceptată în practică ca şi una de compromis, este oferită de condiţia de maximizare a raportului semnal/zgomot la ieşirea filtrului de recepţie 13 în cazul transmiterii doar a unui singur simbol de informaţie Acest lucru presupune: a1 = 0 () () ω=ω⇒GS (4 13) E ≠∀= kk,a00 ce, ţinând cont de ec (4 12), implică: 2 ()()() () ω⋅=ω⇒ω⋅=ωGkGGkG (4 14) EER Condiţie satisfăcută, cu k=1, de alegerea: ()() () ω=ω=ωGGG (4 15) ER Relaţia (4 15) nu impune nici o restricţie deosebită caracteristicii de fază a celor )() ωG respectiv ωG, singura condiţ două funcţii răspuns în frecvenţă, ( ERie ce trebuie să fie verificată este acea de sumă nulă a celor două faze la orice ) ωG reală , cum este şi cazul nostru, putem alege: frecvenţă Pentru o funcţie ( () () () ω=ω=ωGGG (4 16) RE 4 1 Procedura de simulare a lanţului de transmisiuni în banda de bază 4 1 1 Discretizarea semnalului de intrare (at, prin eşantionare cu pasul ) l T Se obţine astfel secvenţa [an anT= Alegând TT L l== >20,,, () e]deeb n rezultă o lungime totală MNL= a secvenţei Alegând ,N2= va rezulta o b⋅ lungime totală a secvenţei, M, tot o putere întreagă a lui 2 4 1 2 Echivalarea filtrelor în timp continuu, dimensionate conform ec (4 16), cu filtre în timp discret Echivalarea se face pe baza metodei de echivalare prin invarianţa răspunsului la impuls Se găsesc funcţiile răspuns în )() ΩG şΩG, frecvenţă ale filtrelor numerice echivalente, ( i, respectiv EdRd l 20, (ce respectă teorema eşantionării), funcţii ce, în condiţiile TT L l== eb > au expresia: Ω L⎧ 2 bπ ⋅ =≤Ω cosL⎪ 0 bΩ () () =Ω=Ω ,LGGG= () 4(4 17) ⎨ bdddRΩ E ⎪ π<Ω<Ω 00 ⎩ () ( )()() π+Ω=Ωπ+Ω=Ω22GG,GG RREEdddd În vederea reprezentării lor pe calculator funcţiile se discretizează, eşantionând 2 ⋅=Ω, rezultând secveţele numerice periodice uniform în punctele kπ Mk 14 [][] ()()GkGΩ= = şi , perioada lor fiind M (o putere întreagă kEddGkGΩkdd ERR a lui 2) Memorarea lor implică reţinerea a doar M valori succesive, (o perioadă completă), din fiecare ] 4 1 3 Se gă seşte secvenţa [ns, de la ieşirea filtrului de emisie Aceasta d se calculează în domeniul Fourier, prin: - găsirea secvenţei [] A k FFT a n=, calculată în M puncte, (de fapt se []() dd 2 ] () obţine secvenţa eşantioanelor [ =, kπ⋅=Ω kddAkAΩk, ale Transformatei M Fourier în Timp Discret a secvenţei numerice [ an A,Ω); () ] dd ][][] - calculul produsului [ kGkAkS⋅=, tot în M puncte; ddd [] () kSIFFTns= - calculul secvenţei [] dd Se obţine astfel semnalul ce urmează a fi transmis prin canal, semnal ce trebuie însumat cu zgomotul canalului 2 ] 4 1 4 Generarea unei secvenţ e de zgomot, [nz, de putere σ, ce trebuie di ] ns, obţinând: însumată cu secvenţa de semnal util [ d [] [] [] nznsn's+= (4 18) ddd ] seşte secvenţa răspuns [nr de la ieşirea filtrului de recepţie 4 1 5 Se gă d Procedura utilizată poate fi una similară celei descrise în §2 3 (prin calcul în domeniu frecvenţă) 4 1 6 Găsirea eşantioanelor [rrn=, cu nkL=⋅, ce corespunde ] kdkkb rt, de la ieşirea filtrului de recepţie eşantionării semnalului în timp continuu,( ) eşte secvenţa din fig 1, în momentele de timp tkT= k În acest scop se înmulţ rn, cu o secvenţă de “eşantionare”, [pn, de forma: []] d nL⎧ b1M pn= (4 19) [] ⎨ in rest0 ⎩ 4 1 7 Luarea deciziilor asupra simbolurilor a transmise pe baza k antion r cu un eşantioane-lor r În acest scop se compară k fiecare valoare de eşk prag d (1/2 în cazul unei secvenţ e unipolare şi 0 în cazul unei secvenţe bipolare), şi se decide simbolul: daca r d≥ 1 k⎧ $ a= k⎨ 0 daca r d< k⎩ 15 5 Partea experimentală 5 1 Utilizând un program de calcul ş tiinţific (Mathcad sau Matlab), vizualizaţi graficul funcţiei răspuns în frecvenţă a sistemului în timp discret 2 cu kπ⋅=ΩM=128 puncte, şi considerând 16=L echivalent, () kdGΩk, în b M 5 2 Aplicând algoritmul TFRI gă siţi şi reprezentaţi grafic funcţia răspuns la impuls, [] ng, al acestui sistem Notaţi poziţiile şi valorile maximelor locale d ale secvenţei Diferenţiind secvenţa, apreciaţi viteza de variaţie a semnalului în jurul trecerilor prin zero Ce se poate constata ? rind procedura descrisă în secţiunea 2 3 generaţi o secvenţă 5 3 Urmă ] aleatoare de date unipolare, [ na, de lungime M=128, versiunea discretizată a d unui semnal a(t) de N=8 simboluri unipolare 5 4 Gă siţi secvenţa [] ce se obţine la ieşirea sistemului numeric dacă la nr d ] intrarea sa se aduce secvenţa [ na În ce domeniu de valori se încadrează d nr ? Justificaţi valorile [] d 5 5 Stabiliţ i momentele de timp optime de eşantionare, n, şi găsiţi k $ secvenţa de biţi obţinuţi după luarea deciziei, { a Ce se poate constata ? } k Explicaţi importanţa stabilirii corecte a acestor momente de timp Ce condiţie impune acest lucru sistemului de transmisiuni implementat ? i 3 3-3 5 pentru cazul unei secvenţe de date bipolare 5 6 Repetaţ teţi numărul de biţi din secvenţa de intrare, fără modificarea lui 5 7 Creş ei de intrare vom constata o decizie incorectă la bL La ce lungime a secvenţ recepţie? 5 8 Utilizând o funcţie de generare de valori aleatoare generaţi o secvenţă de valori aleatoare, [] nz, de lungime M=128 Distribuţia valorilor este d Reprezentaţi secvenţa sumă gaussiană, cu o dispersie de 50 = iσ [] [] [] nznsn's+= ddd 5 9 Repetaţi 5 1-5 7 pentru cazul unui semnal de N=512 simboluri binare aleatore Modificând valoarea dispersiei zgomotului, schiţaţi grafic evoluţia probabilităţii deciziei incorecte în funcţie raportul semnal/zgomot Se vor σ considera valorile 5 1,25 1,1,75 0= i 16